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We describe the effect of hydrodynamic interactions in the sedimentation of a pair of inextensible semiflexible 
filaments under a uniform constant force at low Reynolds numbers. We have analyzed the different regimes and 
the morphology of such polymers in simple geometries, which allow us to highlight the peculiarities of the 
interplay between elastic and hydrodynamic stresses. Cooperative and symmetry breaking effects associated to 
the geometry of the fibers gives rise to characteristic motion which give them distinct properties from rigid and 
elastic filaments. 

PACS numbers: 87. 15. He Dynamics and conformational changes of biomolecules 05. 10. -a Computational methods in statis- 
tical physics and nonlinear dynamics 05.65.+b Self organized systems 87.15.Vv Diffusion of biomolecules 



I. INTRODUCTION 

The understanding of the hydrodynamics of semiflexible 
mesoscopic filaments has gained interest due to the relevance 
of these fibers in different contexts. Many biopolymers are 
virtually inextensible semiflexible and their dynamics in a 
fluid plays a central role in the motion of cilia, and eukary- 
otic and prokaryotic flagella (lh. Although cell motility has 
been investigated for decades [2], recent advances in micro- 
fabrication and micromanipulation enable us to interact di- 
rectly with them in simplified in vitro environments, where 
most of the parameters are under control. This allows direct 
and well defined measurements. For example, Riveline and 
coworkers JH have employed optical tweezers to periodically 
oscillate actin filaments connected to micron-sized beads, in 
order to devise an artificial "one-armed swimmer". More re- 
cently, Dreyfus et aJ fl have been able to produce artificial 
swimmers out of polymer-linked magnetic beads. This ap- 
proach enables an easier control of the filaments through mag- 
netic fields, and has allowed to perform quantitative measure- 
ments of the physical properties of the chains, such as their 
bending stiffness, opening a new method to induce proper- 
ties of the linker molecules J5|] or the affinity of the chemi- 
cal contacts between polymer and particle coating from sim- 
ple mesoscopic measurements, such as image analysis from 
video-microscopy @]. Many other possibilities remain to be 
explored, ranging from the use of semiflexible filaments in mi- 
crofluidic devices to the fabrication of synthetic ciliary arrays, 
or to technological applications of artificial swimmers. 

These advances give a renewed stimulus to a quantitative 
and careful analysis of the hydrodynamics of semiflexible fil- 
aments 0, H [Hi • They differ from flexible polymers in 
the ways in which elastic and hydrodynamic stresses com- 
pete, and it is necessary to treat both on the same footing, 
a theoretical challenge. To gain understanding in this inter- 



play and its dynamic implications it is useful to consider the 
motion of semiflexible chains subject to a uniform external 
driving. A single sedimenting filament has been considered 
recently, and it has been shown theoretically ll U Il2l 1 1 3fl that 
the chain response differs qualitatively from that of a rigid rod, 
in accordance with predictions on collective properties of fiber 
suspensions [14]. Specifically, the inhomogeneous hydrody- 
namic stress along the fiber induced by hydrodynamic inter- 
actions (HI) leads to filament bending and orientation trans- 
verse to the applied field 111 XII . Upon increasing the driving, 
the shape of the filament changes and becomes eventually un- 
stable; the filament then sediments without reaching a steady 
state. 

We will analyze the interactions between a pair of sedi- 
menting filaments, and will study how the combined effect of 
HI and elasticity induce cooperativity in their motion. The re- 
sponse of the fibers depends on the specific geometry; in par- 
ticular the translation-rotation coupling is sensitive to the sym- 
metry of the relative positioning of the fibers. Although it is 
known that translation-rotation couplings lead to an intricate 
behavior in the sedimentation of rigid rod suspensions IU5I1 . 
flexibility leads to new scenarios. 

Hence, we will consider three simplified situations. In or- 
der to address the role of hydrodynamic cooperativity in the 
absence of symmetry breaking and induced rotation, we will 
analyze first two cases where such a coupling is prevented. We 
are then able to find the relevant scaling regimes for the veloc- 
ity and the short- and long-time deformation amplitudes, as a 
function of the inter-filament distance. Subsequently, we will 
focus in two collinear filaments, the simplest geometry where 
rotation is induced. The proposed situations can be realized 
experimentally in a straightforward way, to test our predic- 
tions. 

The remainder of the paper is organized as follows. In sec- 
tion [II] we present our computational model, defining the rel- 
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evant parameters of the system. In section [TTT] we summarize 
the relevant features of single filament sedimentation, which 
will be useful in subsequent sections. In section HVl we ana- 
lyze the sedimentation of a pair of semiflexible filaments. We 
conclude in section [V] with a discussion of the main results 
and their implications. 

H. MODEL 

We study numerically the dynamics of inextensible semi- 
flexible filaments of length L, which are characterized by their 
bending energy. Such an approach is relevant for a large class 
of biological and non-biological polymers including DNA, cy- 
toskeletal filaments and carbon nanotubes 11611 . as well as for 
filaments where the degree of extensibility is negligibly small. 

A semiflexible filament can be described by the arclength 
distance along the filament at a given time t, r(s,t), where 
s e [0,1]. Accordingly, its local curvature is C[r(s, t)] = 
\d 2 r(s, t) /ds 2 \, and once the inextensibility constraint is en- 
forced, the elastic energy is given by the Hamiltonian 

H =l K J C[r(s,t)} 2 d S , (1) 

where k stands for the filament's stiffness. 

We model such a filament as a chain of N spheres ("beads") 
of radius a, connected by bonds of fixed distance b. Corre- 
spondingly, the bending energy is expressed as the discretiza- 
tion of the Hamiltonian, 

N-l 

Ub = lY. (i -«*»*<) . ( 2 ) 

i=2 

where 0i is the angle between the bond that connects bead 
i — 1 to bead % with the one that connects bead i with bead 
i + 1. The need to resolve the conformational change of the 
filament makes the simulations computationally much more 
intensive than when considering rigid rods lfl5ll . The change 
in bending energy Hb due to the change in position of bead 
i determines the bending restoring force acting on it, F,b. 
Inextensibility implies that the total polymer length is fixed, 
L = b(N — 1); this quantity is kept constant by constraint 
forces, Fic, applied at every time step on each bead i I a, 1 1711 . 

It is usual to find biopolymers in suspension. Accordingly, 
we need to account for the interactions with the embeding sol- 
vent. Since the Reynolds numbers are small (in water suspen- 
sions, Re ~ 10~ 6 for micron size filaments moving at mi- 
crometer per second, and Re ~ 1 only for millimeter size 
filaments displacing at millimeter per second) we need to ac- 
count for the coupling between elastic and viscous stresses 
acting on the chains. To this end, we consider that each bead 
i is subject to a local friction force 

F iF =- 7o ( Vi -vf), (3) 

where 70 = 6irr]a is a friction coefficient related to the bead 
size and the solvent viscosity, 77, while stands for the ve- 
locity of the solvent generated by the forces that the rest of the 



beads exert on the fluid at the position of bead i. These dis- 
sipative forces couple hydrodynamically all the chain beads 
through the solvent, giving rise to the hydrodynamic interac- 
tions (HI). According to the standard procedure in polymer 
physics, we describe the flow generated by the filament at the 
level of the Oseen approximation. Hence, the induced veloc- 
ity at the position of bead i can be expressed as lfTTl[l8ll 

where a determines the hydrodynamic coupling relevance 
while e,y = (r^ — ij )/ is the unit vector joining beads i and 
j,with rij = |rj — rj\ being the distance between them; Fj 
refers to the total force acting on bead j. Although alternative 
and more accurate approximations to the induced velocities 
can be implemented [19], this simple coupling is enough to 
capture the essentials of HI on elastic filaments, although for 
small filament separations our prediction will not be in gen- 
eral quantitatively accurate. Following the usual approach in 
polymer physics, we take b = 2a, consistent with the Shish- 
kebab model. The friction coefficient job/ a provides a means 
to relate the model parameters to physical units. The local 
friction force gives rise to an effective filament friction coeffi- 
cient which depends on filament configuration; therefore, this 
approach goes beyond resistive force theory SIM], which re- 
gards the solvent as a passive medium that exerts a constant 
friction coefficient and which does not account for the change 
in friction with filament configuration. 

The description in terms of a local friction force allows to 
describe the filament's dynamics using a molecular dynam- 
ics approach based on the total force acting on each bead, 
F s: = F lB + F iC + F^ + F ie + F ith , with F ie and F ith the 
external force and the random force which accounts for ther- 
mal fluctuations respectively. Multiple filaments can be also 
analyzed without any further algorithmic complexity. In this 
paper we will concentrate on filament sedimentation, where 
Fi e = F e is a constant external force 12111 . in situations where 
the energetic contributions due to the elastic energy domi- 
nate and will therefore neglect thermal forces (Fith = 0). 
Non-uniform oscilatory external drivings have also been con- 
sidered within the same approach in the context of filament 
swimming @, [13, miH. 

Note that the choice b = 2a implies that, as the number 
of beads increases, the hydrodynamic aspect ratio of the fila- 
ment, a/L, decreases as 1/N. The basic hydrodynamic cou- 
pling we analyze depends only marginally on this choice. A 
quantitatively more accurate description of the filament finite 
thickness would require a computationally costly description 
of the filament morphology, and would not affect qualitatively 
the results. We have checked that if we change the number 
of beads , N, but keep the mass and length of the filament 
constant, the results obtained change less than a few percent 
in the worst case for values ranging between N = 30 and 
N = 100. Moreover, we observe a clear tendency to con- 
vergence on increasing N which implies that the trends de- 
scribed subsequently are robust. The results we will describe 
have been obtained for filaments with N = 30. We have 
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taken L — 1 and the filament density p = 1 while the fric- 
tion coefficient 70 = 5 sets the time scale. The bead mass 
is then nib = ir/2N 3 . Numerically, one has to resolve the 
non-physical inertial time scale in which the bead's accelera- 
tion decays; hence we take the time step 1CP 6 . The bending 
rigidity was varied to control the filament flexibility, but the 
characteristic associated time scale is always larger than the 
inertial time scale, ensuring that inertia becomes irrelevant at 
the scales in which the filament configuration evolves. 



III. SEDIMENTATION OF A SINGLE SEMI-FLEXIBLE 
FILAMENT 

We briefly describe the main features of single filament 
sedimentation, which has been explored both analytically and 
computationally Ull [I2I [HI], m order to help with the anal- 
ysis in the coming sections. Filament sedimentation can 
be described in terms of the dimensionless parameter B = 
L 3 F e /n, the ratio of the energy imparted by the external force 
and bending energy 111 111 . We disregard thermal effects, which 
we consider subdominant. 

If an external homogeneous force (F e ) is applied transver- 
sally to the filament axis, its shape reaches a steady state as a 
result of the competition between elastic, constraint, external 
and friction forces. Since the friction force acting on beads 
near the chain's ends is smaller than the local friction in their 
center, filaments bend and align perpendicular to the exter- 
nally applied field. For low to intermediate values of B, chain 
sedimentation can be characterized by the bending amplitude, 
A, defined as the distance between the highest and lowest 
beads along the direction of the applied field. For B < 50, the 
filament's amplitude increases linearly; at B ~ 200 a plateau 
is reached, signaling a saturation of the filament deformation 
in response to the applied field, as depicted in Fig.Q] At even 
larger values of B, metastable shapes with two minima are 
observed Hill . 

If the chain is not aligned perpendicular to the applied 
field, the friction force is not balanced and it generates a net 
torque that will align the filament perpendicular to the exter- 
nal force fill [T3l . For weak forcings, where the degree of 
bending is proportional to B, the torque generated is also pro- 
portional to B. However, the time it takes the polymer to ro- 
tate increases as 1/B for weak forcings, leading to a singular 
behavior for a rigid rod (B = 0), in which the filament keeps 
its initial orientation because it takes an infinite time to rotate 
the filament to its perpendicular orientation, a feature that is 
not captured by resistive force theory. Hence, a single elastic 
polymer reacts to an applied field in a qualitatively different 
way than a rigid one. 



IV. SEDIMENTATION OF A PAIR OF FLAMENTS 

We will now consider pairs of symmetric filaments of 
length L and rigidity k, at a distance d and subject to an ex- 
ternal uniform force field F e . We assume that F e is parallel 
to e 2 and that the polymers lie initially perpendicular to the 




FIG. 1: Bending amplitude of a single and a pair of sedimenting 
filaments, separated a distance d/L — 0.1, as a function of the di- 
mensionless parameter B. 

applied field, along the e^. direction. The details of the coop- 
erativity induced by the hydrodynamic coupling are sensitive 
to the initial configuration. To distinguish between different 
effects induced by hydrodynamics, we will consider three dif- 
ferent geometries, as depicted in Fig. [2] which correspond to 
parallel (geometries I and II), and collinear (geometry III) fil- 
aments. Geometries I (Fig.fSJa) and II (Fig.[2]b) preserve the 
mirror symmetry with respect to the filament's center while 
geometry III (Fig. |2]c) will allow us to explore the effects of 
translation-rotation coupling. We will see that geometry I con- 
serves the symmetries of the one-filament case, and geometry 
II breaks the up- down symmetry. 

The presence of a second chain modifies the friction exerted 
on the first filament. As a result, the filament shape and veloc- 
ity will change as a function of the distance between chains. 
Depending on their initial conditions, the presence of a second 
thread can induce rotation of the sedimenting polymer, break- 
ing the mirror symmetry. We will characterize this translation- 
rotation coupling through the deformation asymetry parame- 
ter, 



where — \xk — x m i n | is the distance along the x-direction 
between the fc-th bead and the lowest bead, as shown in 
Fig. |2]a. The parameter D ranges between [—1, 1], and re- 
flects the transverse asymmetry of the filament ends. For a 
single filament, there is no shape asymmetry and D = 0. 

A. Geometry I. Parallel Filaments Under a Force Transverse 
to the Plane They Define. 

The first geometry under consideration involves two fila- 
ments that are parallel and transverse to the external force 
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FIG. 2: Pairs of semiflexible filaments characterized by their stiffness k and length L that are separated a distance d, when they sediment 
under the action of a uniform external force field, F e , for three different geometries; in the three cases we show the relevant parameters. 
I) Geometry I: Parallel filament sedimenting due to a uniform force transverse to the plane defined by them. II) Geometry II: Sedimenting 
coplanar filaments. Ill) Geometry III: Sedimenting collinear filaments 



(Fig. EJa). Due to the initial configuration, they will sediment 
at the same speed with D = and keeping their initial sep- 
aration d. After a short time interval, in which HI propagate 
and the filaments' inertia decays, they reach their steady state 
sedimentation velocity and deform into shapes analogous to 
the ones described for single polymer sedimentation. 

Since in our model HI propagate instantaneously, the sed- 
imentation velocity of the initially straight filaments after 
one time step will deviate from its free draining value vq = 
F e /Njo; in appendixlAlwe compute this initial velocity. Sub- 
sequently, the filaments will deform until they reach a new 
steady state in which they fall down at a different speed. 
Hence, hydrodynamic cooperativity shows up in the degree 
of filament bending and its dependence on chains separation. 

Bending amplitude. The filament deformation can be 
characterized through the same bending amplitude, A, defined 
for single filament sedimentation. For a given value of B, A 
will now depend on filaments' separation, d. We have found 
that A(d) is consistent with an algebraic decay, 

A(d) ~ cT 3 . (6) 

as displayed in Fig. 0a. The dependence of A on B for a pair 
of filaments separated a distance d does not differ quantita- 
tively from the one observed for a single filament if we com- 
pare the filament's shape with equal values of the parameter 
A. The effect of the second filament can hence be understood 
in terms of an effective bending energy. Making use of Fig. Q] 
it is possible to reproduce the filament's shape by identifying 
B e f f (d) once A has been measured. 

Sedimentation velocity. The presence of a second filament 
leads to an increase of the sedimentation velocity, which de- 
cays algebraically down to small distances, 

v(d) ~ d- 1 (7) 

as shown in the inset of Fig. [3]b. Such a functional depen- 
dence derives from the form of HI at Oseen level [|21]. As a 



result of such coupling, two sedimenting filaments affect each 
other at large distances, and the coupling becomes quantita- 
tively relevant at distances of the order of the filament's size. 

For a given separation d, the velocity change due to hydro- 
dynamic cooperativity decreases with B. More rigid filaments 
have a larger filament section exposed to the flow induced by 
the neighbour filament, leading to a larger relative velocity in- 
crease. 



B. Geometry II. Sedimenting Coplanar Filaments 

We consider next a pair of straight parallel chains separated 
a distance d under the action of a uniform external field copla- 
nar and transverse to the two filaments, as depicted in Fig.|2]b; 
the symmetry of the geometry ensures again D = 0. The 
upper chain bends less than the lower one, A\ < A2, and 
sediments faster because it is subject to a smaller drag due 
to the solvent counterflow. Similar phenomena have been re- 
ported for the sedimentation of other flexible objects, such as 
drops i23il24ll . In Fig. [4] we display the relative sedimentation 
velocity, v r = |i>2 2 | — \ v iz\> as a function of a prescribed in- 
terfilament distance, d. The sedimentation velocity increases 
with B until the filament deformation reaches the plateau de- 
picted in Fig. Q] The relative velocity vanishes at B = and 
for larger values of B increases up to a 10%. Accordingly, the 
sedimentation velocity increases until the plateau regime of 
the filament deformation is reached, where the time scales for 
displacement and bending do not differ significantly. Since v r 
vanishes for B = 0, the relative velocity can change signif- 
icantly upon increasing the filaments' flexibility. Moreover, 
since the relative velocity is always negative the two filaments 
will always approach and will eventually collide. We have 
verified that the filaments sedimentation velocity decays as 
<f -1 , while their relative velocity, v r , decays as d~ 2 (Fig. |4j) 
because the leading contribution of 0(d^ 1 ) cancels out ex- 
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FIG. 3: a) Bending amplitude, A, and b) velocity in the direction of the force field as a function of filaments' separation for Geometry I, at 
B = 1 and B = 200. The asymptotic values of A are defined as A = A(d -> oo), are Ao(B = l)/L = 9.43 x 10~ 4 and A (B = 
200)/L = 3.40 x 10 _1 . Velocities are expressed in units of the sedimentation velocity of an isolated filament, vo(B) = v(d — > oo, B). 




d/L 

FIG. 4: Relative velocity of a pair of initially parallel semifiexible 
filaments sedimenting on top each other (Geometry II) at B = 10 
normalized by the sedimenting velocity of a single filament at the 
corresponding B, vo(B). a) Relative velocity for B — 200, where 
the values of v r /vo(B) are much larger, b) The algebraic depen- 
dence of the relative velocity on the distance is compatible with d~ 2 . 



actly. Such a behavior is general and valid for all values of B 
and in AppendixlBlwe discuss such a dependence on the basis 
of a simplified limiting model. 



C. Geometry III. Collinear Filaments in a Transverse Field 

Finally, we analyze the sedimentation of two collinear fil- 
aments under the action of a uniform transverse field. To 
this end, we consider a pair of filaments which are initially 
straight and with a minimal bead-to-bead distance d, as shown 
in Fig.|2]c. 

The hydrodynamic coupling induces a sedimentation veloc- 
ity which differs from free draining motion, vq. Due to the 
instantaneous propagation of HI in our model, deviatons from 
vq are observed after one time step. In Appendix [? ] we 
compute this initial velocity when bending is negligible. The 
presence of the second filament induces in general a relative 
displacement of the filaments and also a rotation because the 
mirror symmetry is lost. The time scales at which filaments 
rotate and displace depend on filament flexibility. 

At short times, when filaments have deformed significantly 
although their distance has not changed appreciably, one can 
characterize the filaments by a sedimentation velocity, v s , 
which can be understood as the limit v s = lim t ^o v(t) ^ vq. 
At long times, filaments approach or move apart and rotate 
significantly only after they have developed a well-defined 
bent shape. 

In Fig.|5]we show the sedimentation velocity of a pair of fil- 
aments at different values of B, as a function of time in units 
of t c = L£j_/F e , the time it takes a filament to displace its 
own size, where £j_ = Air-qL/ log(L/6) is the friction coeffi- 
cient of a rigid rod in the slender body limit. One can see how 
v s is reached on time scales of order r c , and that a second, 
smaller velocity is reached at larger times, when filaments dis- 
place and B is not too large. Upon increasing B, this second 
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FIG. 5: Sedimentation velocity of a pair of filaments at d/L = 0.5 as 
a function of time. At B = 5 the velocity reaches a plateau at small 
times which is approximately stationary, whereas for B — 150 the 
velocity decays after the initial plateau at t/r c ~ 1. This dependence 
increases with B, as it becomes clear for the curve at B = 1500. 
Inset: velocity increase at short times. 

plateau becomes a slow decay toward the long-time regime. 
The range of this decay also decreases with B. 

In this geometry the long-time sedimentation regime is 
characterized generically by a coupling between translation 
and rotation. We will first describe the short-time sedimenta- 
tion regime, and address the long-time behavior subsequently. 

1. Short times 

The presence of a neighbouring filament induces a defor- 
mation that is not symmetric with respect to the center of mass 
of each filament. The deformation asymmetry increases as the 
filaments approach, with an algebraic dependence 

D(d)~d- 2 , (8) 

as displayed in Fig. |6]a. The change in D arises either be- 
cause of filament tilt, at small B, or due to bending, at large B 
values. Analogous to the observations in the previous geome- 
tries, A(d) decreases algebraically with filament separation, 
as shown in Fig. [6jb, although with an exponent —2 instead 
of —3. The short-time sedimentation velocity v s (d) decays 
algebraically toward the single filament value algebraically as 
1/d. However, the dependence on B is much weaker than 
in previous geometries, as displayed in Fig. |6]c. Thus, in the 
short-time regime we recover practically a universal depen- 
dence of the velocity on the distance. 

2. Long times 

In order to analyze the behavior at long times, it is use- 
ful to analyze separately the small B (B < 200)and large B 



(B > 200) regimes, corresponding to filaments which do not 
reach and reach the saturation of single filament deformation, 
respectively, as shown in Fig.Q]). 

Small B In this regime, two collinear parallel filaments 
always tilt, and rotate as a result of the inhomogeneous hy- 
drodynamic stress along every filament. As a result of the 
rotation and the geometrical asymmetry, we observe a rela- 
tive velocity, v x . In Fig. [7] we display the angular velocity, 
u, which shows a crossover from a d~ 2 dependence on fila- 
ment separation to a weaker \j\fd at shorter distances. On 
the contrary, v x decreases as 1/d. The weaker dependence at 
short distances of the angular velocity develops as a result of 
the translation-rotation coupling; at a fixed distance the angu- 
lar velocity decays always as 1/d 2 . The weaker dependence 
on filament distance and small magnitude of u> implies that 
filaments will usually approach and collide before they have 
rotated by an angle larger than ir/2, which would allow them 
to move away from each other (see the example in Fig. [H]a). 
Only initially remote filaments will avoid collision on observ- 
able time scales. 

Filament rotation can be clearly analyzed if the initial dis- 
tance between the two filaments, d, is fixed. One can observe 
that the two filaments rotate with an averaged angular veloc- 
ity which increases with decreasing d, as seen for example 
in Fig. [8]b. In Appendix ICl we compute the angular velocity 
for a pair of paralell sedimenting filaments, which shows the 
relevance of the translation-rotation coupling. 

Large B When the degree of deformation is limited by 
the length of the filament, its bending is less sensitive to the 
presence of a second fiber. As a result, the angular velocity 
that characterizes rotation decreases significantly with B. As- 
sociated to this reduced sensitivity, we observe that for any de- 
gree of flexibility, at long times there exists a threshold B* (d), 
such as for B < B* the two filaments approach, while for 
B > B* move apart after a transient induced by their initial 
condition, as summarized in the B — d diagram in Fig.[9]a; as 
B increases filaments move apart at smaller distances. Such 
a behavior is not present for rigid filaments, and it correlates 
with the component of the velocity along the filament axis. 
In Fig. |9]b we show such a velocity as a function of B for 
a given initial separation, on a time scale in which filaments 
have displaced distances comparable to their sizes. One can 
see that the velocity reverses sign at a finite value of B, which 
in fact coincides with the change in behavior displayed in the 
B — d diagram. When arriving at the bending plateau, for 
B > 200, filaments move apart (Fig.[l0]a). On the same time 
scales we have computed the degree of asymmetry D, as de- 
picted in Fig. [9J C - The behavior is qualitatively similar to the 
one observed at short times. As shown in the inset, for very 
small values of B the filaments essentially only rotate, and 
only for B > 1 flexibility starts to affect D quantitatively. 
Such parameter D also reverses sign at a value of the flexibil- 
ity similar to that characterizing the change in the velocity. We 
attribute this change in trend to a crossover from a regime of 
small B, where filaments essentially rotate rigidly, to a regime 
where the asymmetric deformation is controlled by bending. 

At large values of B, a third regime is observed both in v x 
and D; both quantities reach a minimum and decrease again 
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FIG. 6: Different measures of the configuration of a pair filaments in the short time regime at B = 1 and B — 200 for Geometry III. In 
all the insets we show the same curves in log-log scale to display their power law dependence which does not depend on B. a) Deformation 
asymmetry, D(d). b) Bending amplitude, A(d), normalized by the corresponding one filament value, A . c) Velocity in the direction of the 
external force. The asymptotic values are vq(B = 1) = 0.52L/r c and vo(B = 200) = 0.50L/r c . 




d/L 

FIG. 7: Angular velocity, u>(d), of a pair of filaments as a function of 
the distance, for different initial separations, at B = 1 for Geometry 
III. u){d) is normalized by u>o, the angular velocity at the closest dis- 
tance. The different curves correspond to different initial separations 
(d(t = 0) = 5, 3, 2, 1 starting from right to left); after a transient 
related to the relaxation of the (initially straight) filaments, they fol- 
low a unique curve which depends only on the distance betwen the 
filaments. The inset shows the algebraic behavior at small and large 
distances. 



in magnitude. The minimum is observed in the parameter re- 
gion where metastable filament configurations for single fila- 
ment sedimentation develop. At even larger values of B, we 
have also observed regimes where intrafilament collisions are 
observed as a result of the shape deformations the filament 
suffers during its sedimentation. Therefore, the final configu- 
ration is not stable and changes continuously with time. 



V. CONCLUSIONS 

We have studied the sedimentation of a pair of filaments 
suspended in a low Reynolds number fluid. The coupling 
that the filaments in the solvent induce on each other through 
flows, the so called hydrodynamic interactions (HI), give rise 
to a rich variety of cooperative motion. We have concentrated 
on the simplest geometries, in order to perform a careful anal- 
ysis that allow us to focus on the essential features of such 
cooperativity. To this end, we have implemented and used a 
simple and efficient numerical method which models the fil- 
ament as a set of beads and imposes inextensibility where HI 
are treated at the Oseen approximation. 

The geometries we considered have helped us to show that 
in all cases the sedimentation of a pair of semiflexible fila- 
ments is qualitatively different from that of rigid filaments. 
The rigid limit is in fact singular, and sets in because the 
time the filaments need to modify their initial configuration in- 
creases with filament rigidity, and diverges for infinitely rigid 
rods. For sufficiently symmetric geometries, such as Geom- 
metry I, HI modify the degree of deformation of each filament 
and its final sedimentation velocity. The interaction decays 
algebraically with filament distance, and it becomes quanti- 
tatively relevant for separations of the order of the filaments' 
size. 

For sedimenting parallel coplanar filaments, we have shown 
that the top filament bends more and moves faster, inducing 
pair collision, as opposed to the sedimentation of rigid rods. 

In less symetric situations, the hydrodynamic coupling in- 
duces a rotation and translation of the sedimenting filaments, 
because their bending lacks the symmetry with respect to the 
filaments' center of mass. In these cases, pairs of filaments 
still interact at long distances, but their sedimentation behav- 
ior becomes more involved, and depends on their degree of 
flexibility as well as their initial conditions. We have shown 
that such rich behavior includes periodic bound trajectories, 
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FIG. 8: a) Configurations of two initially straight filaments a distance d/L = 1 away from each other and B = 100 shown at intervals 
t/T c = 1.2. At long times (t/r c >> 50) the two filaments eventually collide, b) Configurations of two initially straight filaments a distance 
d/L — 0.4 away from each other and B — 1 shown at intervals t/r c — 1.1. In this case we have fixed d to highlight the rotation of the 
filaments, characterized by a mean angular velocity which depends on the relative configuration of the two filaments. The apparent shortening 
of the filament is due to the different scales in the two axis where distances are expressed in units of the filament length L. 




FIG. 9: Effective attraction and repulsion in geometry III. a) B — d diagram showing the distinct behaviors of sedimenting collinear filaments 
as a function of initial distance and driving strength, b) Relative velocity at t = 20r c , when d = L, as a function of B; positive velocity means 
attraction. Inset: decay of the relative velocity at small values of B. c) Deformation asymmetry at d = L as a function of B, showing the 
change of sign and the minimum value at B ~ 3000. Inset: decay at small values of B. 



filament rotation as well as sedimentation with unsteady con- 
formations. 

To sum up, filament flexibility and hidrodynamic coupling 
modify profoundly the behavior of filament sedimentation; 
the simplified geometries explored have helped to understand 
the interplay between elasticity and hydrodynamics and opens 
the possibilties to analyze in detail how such interactions mod- 
ify the response of filament suspensions to applied external 



fields. The study we have carried out has allowed us to make 
definite predictions that can be tested in controlled experi- 
ments, for example with a centrifuge coupled to optical mi- 
croscopy. 
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FIG. 10: a) Configurations of two initially straight filaments a distance d/L = 0.2 away from each other and B — 250 shown at intervals 
t/r c = 2.0. Filaments repel each other and move apart with a highly curved shape, as shown in the inset for the left fiber, b) Configurations 
of two initially straight filaments a distance d/L — 0.2 away from each other and B = 10 4 shown at intervals t/r c — 2.3+. Filaments 
initially evolve into the metastable configuration described in Fig. Q] while they move apart from each other. During that separation the 
filaments change their shape and recover transiently this metastable shape. In the inset we show one of the transient, deformed configurations 
the filaments explore. 



Acknowledgments 

We are grateful to H.A. Stone, for pointing out to us the 
interest of geometry II. I. LI. and LP. acknowledge financial 
support from DGCIYT of the Spanish Government (FIS2005- 
01299). LP. thanks Distincio de la Generalitat de Catalunya 
for financial support. 



APPENDIX A: INITIAL SEDIMENTATION VELOCITIES 

In the initial stages of their sedimentation, straight filaments 
have not deformed significantly. In this regime it is possible 
to obtain analytical expressions for their sedimentation veloc- 
ities at the Oseen level. 

In particular, we are interested in the velocity that a straight 
filament oriented along the x direction induces in a second 
collinear filament a distance d away, as depicted in Fig. [2]c. 
The external field is applied perpendicular to both filaments, 
along the z direction, and hence the distance between beads 
reduces to their separation along the x direction. As a result, 
the velocity on bead i due to bead j in the direction of the 
external force, can be expressed as 



3a F e 
A 70 



i 



(Al) 



where we have assumed that bead i belongs to the filament on 



the left while j is a filament belonging to the filament on the 
right hand side of the pair, and hence the sum runs over the 
beads of this second filament. If we approximate the sum by 
an integral over the length of the filament, we arrive at 



,H 



3a F e 
AL l0 



2L+d 



dx 



L+d 



f _ 3a F e 1 (2L + d 
Xi\ AL 70 \ L + d — x 



(A2) 

The contribution of the second filament to the sedimenta- 
tion velocity of the first one is obtained by computing the 
induced center of mass velocity. In the continuum approxi- 
mation, this induced velocity can be written down as = 

1/L J vj[l dxi, leading to 



3a F e 
AL 70 



In 



d{d + 2L) 



(d + L) 2 



2 In 



d + 2L 



d + L 



(A3) 

We can proceed analogously, with obvious modifications in 
the geometry, to obtain the contribution of a second filament to 
the sedimentation velocity of the reference one for two copla- 
nar filaments (geometry II). In this case we arrive at 



,H 



3a Fe I / d 2 + ^ 2 
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-L + Vd 2 + L 2 
(A4) 

Finally, for Geometry I, if we take the axis of the filaments 
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FIG. 11: Configuration of a pair of trimers under the action of an 
external force. The difference on bending amplitude is given by the 
parameter e. 



along the x axis and their distance along the z we get, 



,H 



- J 2 — 

~ 4L 70 I L 



ct 2 + L 2 



L 2 



ln 




In geometries I and II the velocity diverges as the two fila- 
ments approach each other because all the beads become in- 
finitely close, as opposed to geometry III which is charac- 
terized by a finite induced sedimentation velocity at contact, 

31n(V2)aF e /i7o. 



APPENDIX B: TRIMERS IN GEOMETRY II 

In Appendix[A] we have computed the initial sedimentation 
velocity for a pair of filaments where one of them moves on 
the wake of the neighbouring one (eq. lA51 l. Here we provide 
an estimate of the rate at which they approach each other when 
filament deformation is small (either because we focus at short 
times, or because B is small). We consider the simplest pos- 
sible case, where the filaments are represented by trimers. 

We will calculate the hydrodynamic velocity due to the 
beads of the neighbor filament on the central beads (depicted 
in black in FigfTTTi. defined as c\ and C2 in Fig[TT| The differ- 
ence between the two velocities is a measure of their relative 
velocity. 

We characterize d as the distance between the corner beads 
and the bending amplitude as the separation between the cen- 
tral and corner beads of each filament along the direction of 
the applied field. The top and bottom chains will bend an 
amount A\ and respectively, where = A\ + e. Hence, 
e measures the relative bending amplitude due to the different 
hydrodynamic coupling. The distance between consecutive 
beads of a given trimer is 6 = L/2, and the force field has a 



magnitude F e . 

The velocity on the central beads 1 and 2 has only compo- 
nent in the z direction, 



2{d-A 1 f + b 2 -{A 1 +e) 7 
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"2" To" \[(d-^ 1 ) 2 +6 2 -(A 1 +e) 2 ] 3/2 ' d + ej 
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2 7o { [(d+ Ai +e) 2 + 6 2 - A\f' 2 d + ej' 

(B2) 

Therefore, the relative velocity, v r = v c % — v C 2, can be ex- 
pressed as, 

3a F e f 2(d - Ai) 2 + b 2 - {A l + e) 2 
2 7o \ [(d - Ax) 2 + b 2 - (A x + e) 2 ] 3/2 
2(d + A 1 +e) 2 + b 2 -A 2 \ 



[(d + Ax + e) 2 + b 2 - A 2 } 



2l3/2 



(B3) 



For small B, or large distances, also the differential bending 
will be small. In this case, if we expand the previous expres- 
sion in powers of e, we arrive at 



v r (d, e) 



3aF e 
7o 



(2A 1 + e + e 2 )i_ + 0(e 3 ,d- 4 ), (B4) 



which we have validated using simulations of trimers in this 
geometry. This expression shows that the top trimer moves 
faster as a result of HI. The relative velocity depends on the 
distance as 1/d 2 . The increase of v r with Ax indicates that 
also the relative velocity will increase with the flexibility B. 



APPENDIX C: ROTATION OF FILAMENTS 



In section lTVC2l we have seen that in geometry EH, small 
B filaments rotate and that the parallel filament geometry is 
unstable. In this appendix we consider a pair of straight fila- 
ments, with vanishing bending and constraint forces, oriented 
along the direction of the external field, F e , separated a dis- 
tance d on the x direction. Hence, we can compute the veloc- 
ity on each bead due to the presence of the second filament 
following the same approach as in Appendix [A] The com- 
ponents of the velocity on monomer i of a given filament in 
the directions along and perpendicular to the filaments can be 
expressed as 



(-a = 



3ci F e Xij Zij 
46 70 I-*. 

3a F e x ^ 
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1 + (zij/nj) 



i/b 



(CI) 



where x* 



x, 



d , = \zi — Zj\ and 77 



x 2 j + z 2 - , and where the sums run over the monomers of the 
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same and the neighbouring filaments. Due to the symmetry 
of the configuration, {v^) x only has nonvanishing contribu- 
tions from the neighbouring filament, while (vf) z has a con- 
tribution for the filament itself, which corresponds to the sed- 
imentation velocity of an isolated filament aligned parallel to 
the external field, and which in the continuum approximation 
leads to a sedimentation velocity (v^)f — F e \n(L/b)/2-KriL 
in the slender body limit. The contribution of the neighbour- 
ing filament increases this velocity to, 



verse velocity. Due to the symmetry, this velocity dis- 
tribution does not lead to any transverse translation; in 
fact, (v?) x = (l/£)Ej(*>f)s = 0. This veloc- 
ity profile yields a net rotation of the filament, which is 
given, at large distances by Ioj = m^,^ Zij(vf) x = 

(3amF e /4- f0 L) [Z/d - 7i 2 /12d 4 ] , where I = mL 2 /12 is 
the moment of inertia, hence the angular velocity decays 
asymptotically as 1/d 2 , i.e. faster than the approaching ve- 
locity (uf ) z , 



- + - H l + V^TT 2 ) 



(C2) 

While the filament sediments, it will experience a trans- 
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